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We consider the problem of superconducting Ginzburg-Landau (G-L) vortices with antiferromag- 
netic cores which arise in Zhang's SO(5) model of antiferromagnetism (AF) and high temperature 
superconductivity (SC). This problem was previously considered by Arovas et al. who constructed 
approximate "variational" solutions, in the large n limit, to estimate the domain of stability of such 
vortices in the temperature-chemical potential phase diagram. By solving the G-L equations numer- 
ically for general k, we show that the amplitude of the antiferromagnetic component at the vortex 
core decreases to zero continuously at a critical value of the AF-SC anisotropy (g w 0.25) which is 
essentially independent of k for large k. The magnetic field profile, the vortex line energy and the 
value of the B-field at the center of the vortex core, as functions of anisotropy are also presented. 
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I. INTRODUCTION 

Zhang 10 has proposed a model which attempts to 
relate the antiferromagnetism and d-wave superconduc- 
tivity observed in the high T c cuprates by a symme- 
try principle. In this model, the Neel vector, m, which 
is the order parameter for antiferromagnetism, and the 
real and imaginary parts of the complex superconduct- 
ing order parameter vb form the five components of a 
so-called "superspin" vector which transforms under the 
group SO (5). This symmetry is explicitly broken by an 
anisotropy which favors antiferromagnetism (AF) for the 
undoped compounds (which have one unpaired electron 
spin per site). The anisotropy can be varied by a chemical 
potential term which adds holes (or equivalently removes 
electron spins). For sufficiently large chemical potential, 
the anisotropy changes sign stabilizing the superconduct- 
ing (SC) state. 

A straightforward consequence of this model, as noted 
by Zhang Q, is that, close to the SC-AF phase bound- 
ary, vortices in the superconducting state should have 
antiferrimagnetic cores. The idea is simply that, for 
small anisotropy, it is easier to simply rotate the order 
parameter in the core into the AF direction than it is 
to reduce its amplitude to zero. Arovas et al. fol- 
lowing Zhang's suggestion, calculated the approximate 
domain of stability of vortices with AF cores, within 
Ginzburg-Landau theory, in the temperature-chemical 
potential plane. Their calculation compared the ener- 
gies of two variational solutions for the behavior of the 
superspin vector, (1) the function tanh(r/£) for the mag- 
nitude of the SC order parameter, representing a vortex 
with a normal core and (2) a vector making an angle 
9(r) — (ir/2)exp(—r/£) toward the AF direction from 
the SC direction. The contribution of the vector poten- 
tial to the free energy was! ignored, as is appropriate 
for the extreme type-II limit (k — > oo). By comparing 



the energies of these two configurations, minimized with 
respect to £ and £, as a function of the G-L parameters, 
Arovas and coworkers located a boundary within the SC 
region of the phase diagram, where vortex core solutions 
(1) and (2) exchange stability. If these solutions were 
exact, the implication would be that the amplitude of 
antiferromagnetism in the core changes discontinuously 
across this boundary. 

In this paper we consider the full numerical solution of 
the coupled non-linear Ginzburg-Landau equations for 
vortices in the phenomological SO (5) model, including 
the vector potential. We calculate the AF and SC or- 
der parameter profiles and the magnetic field profiles as 
functions of distance from the core for general k and 
anisotropy. We find a much broader range of stability 
for vortices with AF cores than was found by Arovas et 
al. We also find that the amplitude of the AF compo- 
nent in the core falls continuously to zero with increas- 
ing chemical potential as the phase boundary for normal 
cores is approached, and we calculate the critical value of 
the anisotropy(chemical potential) where this transition 
occurs for a range of k, including n — oo. We also obtain 
the dependence of the vortex line energy and the field 
at the core on k in the large k limit for AF and normal 
cores. Our results allow us to conclude that both H c i a! 
nd H C 2 increase with increasing anisotropy as the cores 
become less antiferromagnetic. 



II. FREE ENERGY 

We consider the following Ginzburg-Landau free en- 
ergy for superspin vector n = (tpi, mi, m%, 1^2) with 
vector potential A, 
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The anisotropy is measured by g = 4x(/i 2 — /x 2 ), where /i denotes the chemical potential. 

We introduce a non-dimensional form of the energy. In these units the penetration depth A = 1, and the correlation 
length £ and GL parameter k are related via £ = k^ 1 . More precisely, we let 
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Then the free energy becomes (up to a factor of 2 m*b ) 
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The anisotropy term is regulated by a new constant g 7 
9 4 X (/i 2 - H 2 C ) . . 

.9 = - = n ' ( 5 ) 
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which is positive when superconductivity is preferred in 
the bulk (/i > /i c ). For convenience, we drop the hats on 
non-dimensional quantities in the following. 

A. Vortex solutions 

We now seek vortex solutions in the plane, in other 
words free energy minimizcrs of the form ip(x) = f(r)e l6 , 
m(x) = m(r)m with rh a fixed unit vector, and A(a;) = 
[S(r)/r 2 ](—y, x) = [S(r)/r]6. In addition we require that 
/(0) = 0, f(r -> oo) -> 1, m(r -> oo) -► 0, 5(0) = 
0, and 5(r — > oo) — > 1. Note that requiring m to lie 
in a fixed direction does not restrict the minimization 
problem, since spatial variations in direction of m will 
increase the free energy. With this ansatz, the energy 
takes the form: 
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/, S, m satisfy the following Euler-Lagrange equations: 
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By linearization at r — we obtain the standard prop- 
erties for vortex solutions: HD f(r) ~ ar and 5(r) ~ /3r 2 
near r = 0. Unconstrained minimization also yields 
to'(0) = 0. Using Theorem 7.2 in Jaffe-Taubes |g we 
see that /, S, m approach their asymptotic limits r — * oo 
exponentially. In addition, f(r) and S(r) are monotoni- 
cally increasing, while the ^-component of the local field 
h(r) — S'(r)/r is monotonically decreasing. The AF or- 
der parameter m(r) is either positive and monotonically 
decreasing or identically zero. || These properties are 
illustrated in Figure |l|. 

The results shown in Figure [l] and elsewhere in this 
paper come from numerical solutions of Our 
technique is to find free energy minimizers by computing 
a negative gradient flow for a discretized free energy. We 
approximate the improper integral in E by an integral 
over a finite interval, large compared to the fundamen- 
tal scales (A = 1 and £ = in the problem. We set 
a fixed boundary condition for /, S, m at the right end- 
point, equal to their asymptotic values. We then use fi- 
nite elements to discretize the integral, with a fine mesh 
concentrated at the vortex core. For example, in Fig- 
ure [| we divide the r-interval [0, 10] by 350 grid points, 
with 300 points lying in [0, 2]. We choose standard linear 
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elements, except in the first interval near 0, where the 
behavior of S(r) is known to be quadratic. 

We then explicitly calculate the gradient of the 
discretized energy with respect to the coefficients 
(/fe), (Sk), (wfc) of the base elements. To find the minimal 
configuration, we start with some initial data, and flow 
the coefficient vector (/&), (Sk), (rrik) against the gradient 
until an equilibrium is found. To solve the gradient flow 
numerically, we use a code developed by Dr. N. Carlson 
of Los Alamos National Laboratory, which he very kindly 
made available to us. In this code the ODE is solved by 
an implicit method, using second order backwards differ- 
ences. At each time step the resulting nonlinear system 
is solved by a modified Newton method, and the time 
step is chosen based on estimates of the local truncation 
error. 
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Radial coordinate r Radial coordinate r 

FIG. 1. Typical profiles for f,h,ma,tK— 15. On the left, 
the core is AF, while the core is normal on the right. 



An important difference between our results and those 
of Arovas et al. is that our solutions have superspin mag- 
nitude strictly less than one in the core which allows the 
transition to be continuous. To see this, let u = f 2 + m 2 , 
so u solves the ODE 

u" + ~v! + 2k 2 (1 - u)u = 
r 

2 ^7#^/ 2 + 2 ^ 2 ™ 2 + (/') 2 + K) 2 - (io) 

First, assume that the maximum value of u is achieved 
at tq > 0. Note that u'(ro) — 0, which follows from 
standard calculus if ro > 0, and u'(0) < if ro = 0. In 
addition, u"{tq) < and the right-hand side of ( pp| ) must 
be strictly positive at ro, so 2k 2 (1 — u)u > 0. Since u > 

i 2 < 1. In the 



by its definition, we must have u = j -f 
core region, we can say even more. Since to is positive 



and decreasing, it attains its maximum at r = 0, and 
therefore (||) implies 

> m"(0) = y (g - 1 + to(0) 2 )to(0). 

When m(0) > we have the upper bound on the AF 
order parameter at the core: 

m(r) 2 < m(0) 2 < 1 - g. (11) 

Since f(r) ~ r for r ~ 0, we conclude that the superspin 
magnitude is bounded by ^/l — g at r = 0. 

B. Phase transition 

Conventional superconducting vortices with normal 
cores are obtained by minimizing £ under the constraint 
m(r) = 0. Note that when m = equations re- 
duce to the well-known vortex equations of Abrikosov 
which have been extensively studied in the physical 
and mathematical literature. We denote by 

/k,0j <Sk,o these normal-core vortex solutions. 

We find that there is a critical value g* with 0.2 < g* < 
0.3 at which the stability of the normal-core solutions 
/k,o, £^,0 (with to = 0) changes: normal cores are stable 
(i.e., energy minimizing) when g > g* and unstable for 
< g < g*. To see this, we linearize (||) at the normal 
core solution to obtain an eigenvalue problem: 

J_ m "( r ) + -i-m'(r) + (1 - fl )rh = grh. (12) 

Define g* to be the ground-state eigenvalue, with associ- 
ated eigenstate rh K . From Rayleigh's principle, to k is the 
minimizcr of 
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When g crosses below the value <?*, the linearized oper- 
ator associated with (||) possesses a negative eigenvalue, 
and the normal core solutions become unstable. In fact, 
by perturbing the normal core solution m = by a small 
multiple em of the eigenstate, we see that the energy 
is decreased by allowing a nonzero AF component when 

9 < 9k- 
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Our numerical computations show that the transition 
is continuous (second order) at g*. Indeed, in Figure ^ 
we plot a bifurcation curve, g vs. m(0) obtained by nu- 
merical simulation at several values of n. Since m(r) is a 
decreasing, positive function, the value m(0) represents 
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the maximum value of m. Our analytic and numerical 
results predict that the superspin parameter will lift con- 
tinuously out of the superconducting plane as g is de- 
creased through the critical value <?*, with its maximum 
value approaching 1 as g — > from above. In Figure || 
we illustrate the continuous growth of the AF compo- 
nent with decreasing g. We see also that the radius of 
the vortex core increases as anisotropy g decreases. 



obtained with slowly varying g were checked against solu- 
tions calculated at some fixed g and found to be virtually 
identical. As a further check on our routines we made an 
independent computation of the bifurcation point g* us- 
ing the Rayleigh formula ([l3]). Using a finite element dis- 
cretization of ( |i~3| ) and a standard eigenvalue solver from 
EISPACK, we obtain good agreement with the observed 
transition values in Figure ^. Some values are given in 
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FIG. 2. Nucleation of antiferromagnetic cores for various k 
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20 


0.2502 


40 


0.2531 


80 


0.2541 


120 


0.2543 


oo 


0.2545 



TABLE I. Values of g*. 
1501 grid points. 



Computations were made with 



K=20, g.0. 26, 0.23,0.20,0.1 5,0.1 ,0.05 K.20, g=0.26,0. 23, 0.20, 0.15,0.1 ,0.05 
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Radial coordinate r Radial coordinate r 

FIG. 3. /, m for different values of g. As g decreases, we see 
that the radius of the vortex increases and that m increases. 

In computing the bifurcation curves shown in Figure ^ 
we allow g to change very slowly in time, g(t) ~ 10 -8 t. 
In this way, the minimum point of the free energy moves, 
but very slowly compared to the relaxation time of the 
gradient flow. We ran the program with many different 
values of k and grid spacings and always observed the 
same qualitative behavior for the transition. The results 
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Next we consider the dependence on g of the dimen- 
sionless vortex-line energy ej of an AF core for large k 
(k = 200). A standard estimate (|] shows that for large 

K 



we obtain the extreme Type II equations for 
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Our numerics predict a value C\(g) = 0.49739 when g > 
g^ (normal cores), which is in good agreement with the 
value C\ = 0.4976 obtained by Hu [§. For AF cores, we 
observe that C\(g) decreases as g decreases (c/. Figure ||). 




0.05 0.1 0.15 0.2 
Anisotropy g 



0.05 0.1 0.15 0.2 0.25 
Anisolropy g 



FIG. 4. Dependance of Ci and Co on g 



Similarly, for large k 



h(0)~ In k+C ( 9 ) 



(16) 



with Co(g) = —.28192 for g > g^ (in good agreement 
with Hu || who obtained C = -.2823) and with C (g) 
decreasing as g decreases (cf. Figure ||). So the values 
of ei and h(0) decrease as the cores becomes more AF. 
Recall that the size of the core also increases as can be 
seen in Figure ^. These results show that both H Cl and 
H C2 decrease as the core becomes larger and more anti- 
ferromagnetic. 
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For this model, it is rigorously proven [0 that a con- 
tinuous phase transition to AF cores occurs at a critical 
value g — g*^ , with g* the ground-state eigenvalue of the 
linearized equation (|18|) at the normal core solution, 
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At each value of g there is exactly one stable vortex solu- 
tion. That stable profile has a normal core when g > g^ , 
and an AF core when g < g^. In other words, the bi- 
furcation diagrams obtained by numerical simulation of 
with k finite are proven valid for the limiting 
problem (|l7|), (|l8|). This provides another check on our 
simulations. For k large we should see qualitatively sim- 
ilar behavior between this model and that of ([l7]), (|T§|). 
Numerical computation of g^ using the Rayleigh formula 
dl9| ) is cons isten t with the other calculated values of g*. 
(See Table [IB.) It also shows that vortices with AF 
cores are stable in a much broader range than was found 
in Arovas et al. [^J, where g^ is estimated as 0.0941. 
This discrepancy is due to the fact that the superspin 
magnitude is in fact far from 1 in the core region (see 
our analysis in Section 2.1), while the variational solu- 
tion of Arovas et al assumes that AF core solutions 
have superspins constrained to the unit sphere. 



IV. PHASE DIAGRAM 

The results of the previous two sections for g* define 
a region in the phase diagram in which vortices with an- 
tiferromagnetic cores are stable. To illustrate this state- 
ment, we add the (dashed) line for the stability bound- 
ary of vortices with antiferromagnetic cores to the SO (5) 
mean-field phase diagram, assuming the relation (|^) be- 
tween g and /j,. The result is shown in Figure ^. 



III. THE EXTREME TYPE II LIMIT 



To compare our results more directly with the work of 
Arovas et al Q we consider the extreme type II model. 
In fact, for this simpler model we may obtain rigorous 
analytical results concerning the transition to AF cores. 
H After rescaling by the correlation length £ = n -1 , the 
order parameters converge (uniformly) to a fixed profile, 
/k(^) -* JooM and m re (£) -> ra^r) as n -> oo. (Note 
that in this scaling S K (^) — > 0.) Passing to the limit in 
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Phase diagram with d-0.1 , j-0.9, g -0.2545 



(1.904,0.8227) 

\ 




FIG. 5. Mean field phase diagram for the "spin-flop" model 
defined by Eq. (po|). The solid lines separating the AF and SC 
bulk phases from the high temperature disordered state repre- 
sent continuous transitions. The dot-dash line separating AF 
bulk from "SC bulk AF cores" is first order, and the dashed 
boundary between "AF cores" and "normal cores" represents 
the continuous transition obtained in this paper. 

The SO (5) mean- field phase diagram of Figure || 
is actually the mean field "spin-flop" phase diagram 
for the three-component antiferromagnet with uniaxial 
anisotropy in a parallel magnetic field. In this analog sys- 
tem jy, the antiferromagnetic state with its Neel vector 
parallel to the field corresponds to the SO (5) AF state, 
while the "flopped" antiferromagnetic state, with Neel 
vector perpendicular to the applied field, corresponds to 
the SC state. The field H plays the role of the SO(5) 
chemical potential /i. The spin Hamiltonian for this sys- 
tem is: 
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Where d = D / J and j are parameters which control the 
anisotropy. Figure shows the mean field phase dia- 
gram for this model for specific values of these parame- 
ters (d=0.1, j=0.9). The unit of energy in this diagram 
is To = zJ/3 where z is the coordination number. The 
region labeled "SC bulk, AF cores" is bounded on the left 
by the first order SC-to-AF transition and on the right 
by the stability boundary for AF cores. 

At the mean field level one might expect to see an- 
tiferromagnetic cores throughout this region, with the 
amplitude of the antiferromagnetism being largest close 
to the first order boundary. However, as Bruus et al fl~l| ] 
discuss, the spins in the cores are likely to be only weakly 
correlated from one vortex to the next, and hence the sig- 
natures of antiferromagnetic vortex cores may be rather 
subtle. Bruus et al suggest that excitations of the cores 



could be studied by inelastic neutron scattering, while 
Arovas et al j|] proposed studying the effect of antiferro- 
magnetic cores on the /iSR spectrum. 



V. ACKNOWLEDGEMENTS 

This work was supported by the Natural Sciences and 
Engineering Research Council of Canada. We thank the 
Brockhouse Institute for Materials Research for support- 
ing a workshop which brought together physicists and 
mathematicians to discuss issues in superconductivity. 
AJB acknowledges the hospitality of the Aspen Cen- 
ter for Physics, where part of this paper was written, 
and support from the Superconductivity Program of the 
Canadian Institute for Advanced Research. 



[1] 
[2] 

[3] 

[5] 

[6] 

[7] 

[8] 
[9] 

[10] 



[11] 



Rev. 



S.C. Zhang, Science, 275:1089, 1997. 
D. Arovas, J.Berlinsky, C.Kallin, S.C.Zhang, Phys. 
Lett., 79:2871, 1997 

M.S. Berger, Y.Y. Chen, Jour. Functional Analysis, vol. 
82, 259-295, 1989. 

B. Plohr, Princeton University PhD Dissertation, 1980. 
A. Jaffe, C. Taubes, "Vortices and Monopoles," 
Birkhauser, Boston 1980. 

S. Alama, L. Bronsard, T. Giorgi, preprint 



math.AP/99031271 v2. 

Zh. Eksperim. % Teor. Fiz. 32, 1442 
JETP 5,1174 (1957)]. 

B, vol. 6, 1756-1760, 1972. 



A. A. Abrikosov, 
(1957) [Sov. Phys. 
C.-R. Hu, Phys. Rev. 
M. Tinkham, "Introduction to Superconductivity," 
McGraw-Hill, New York, 1996. 

Arovas et al. sketched a calculation similar to the above, 
for a free energy with anisotropic fourth order terms, in 
Eqs. (8-10) of their paper. However they did not present 
any numerical results. 

H. Br uus, K.A. Eriksen, M. Hallundbaek and P. Hede- 
gard, |cond-mat/9807167|/. 



G 



